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Abstract 

We describe the new version (vl.75r) of the code hfodd which solves the nuclear 
Skyrme-Hartree-Fock problem by using the Cartesian deformed harmonic-oscillator basis. 
Three minor errors that went undetected in the previous version have been corrected. The 
new version contains an interface to the LAPACK subroutine zhpev. Several methods of 
terminating the Hartree-Fock iteration procedure, and an algorithm that allows to follow 
^ ■ the diabatic configurations, have been implemented. 
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§ ■ NEW VERSION PROGRAM SUMMARY 



Title of the program: hfodd (vl.75r) 
Catalogue number: 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ire- 
land (see application form in this issue) 

H 

Reference in CPC for earlier version of program: J. Dobaczewski and J. Dudek, Comput. 
Phys. Commun. 102 (1997) 183 (vl.60r). 

Catalogue number of previous version: ADFL 

Licensing provisions: none 

Does the new version supersede the previous one: yes 

Computers on which the program has been tested: CRAY C-90, SG Power Challenge L, 
IBM RS/6000, Pentium-II, Athlon 

Operating systems: UNIX, UNICOS, IRIX, AIX, LINUX 



1 E-mail : j acek . dobaczewski@fuw . edu . pi 
2 E-mail: j erzy. dudek@ires.in2p3.fr 
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Programming language used: FORTRAN- 77 

Memory required to execute with typical data: 10 Mwords 

No. of bits in a word: The code is written in single-precision for the use on a 64-bit 
processor. The compiler option -r8 or +autodblpad (or equivalent) has to be used to 
promote all real and complex single-precision floating-point items to double precision 
when the code is used on a 32-bit machine. 

Has the code been vectorised?: Yes 

No. of lines in distributed program: 23 987 (of which 10 445 are comments and separators) 

Keywords: Hartree-Fock; Skyrme interaction; Self-consistent mean-field; Nuclear many- 
body problem; Super deformation; Quadrupole deformation; Octupole deformation; Pair- 
ing; Nuclear radii; Single-particle spectra; Nuclear rotation; High-spin states; Moments 
of inertia; Level crossings; Harmonic oscillator; Coulomb field; Point symmetries 

Nature of physical problem 

The nuclear mean-field and an analysis of its symmetries in realistic cases are the main 
ingredients of a description of nuclear states. Within the Local Density Approximation, 
or for a zero-range velocity- dependent Skyrme interaction, the nuclear mean-field is local 
and velocity dependent. The locality allows for an effective and fast solution of the 
self-consistent Hartree-Fock equations, even for heavy nuclei, and for various nucleonic 
(n-particle n-hole) configurations, deformations, excitation energies, or angular momenta. 

Method of solution 

The program uses the Cartesian harmonic oscillator basis to expand single-particle wave 
functions of neutrons and protons interacting by means of the Skyrme effective interac- 
tion. The expansion coefficients are determined by the iterative diagonalization of the 
mean field Hamiltonians or Routhians which depend nonlinearly on the local neutron 
and proton densities. Suitable constraints are used to obtain states corresponding to a 
given configuration, deformation or angular momentum. The method of solution has been 
presented in: J. Dobaczewski and J. Dudek, Comput. Phys. Commun. 102 (1997) 166. 

Summary of revisions 

1. An error in the calculation of one of the time-odd mean- field potentials has been 
corrected. 

2. A factor in the calculation of the multipole moment Q22 has been corrected. 

3. Scaling of the coupling constants has been corrected. 

4. An interface to the LAPACK subroutine zhpev has been created. 

5. Several methods of terminating the Hartree-Fock iteration procedure have been 
implemented. 

6. An algorithm that allows to follow the diabatic configurations has been implemented. 

7. Saving of auxiliary data for a faster calculation of the Coulomb potential has been 
implemented. 
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8. Calculation of average quadrupole moments and radii of single-particle states has 
been added. 

9. Calculation of the Bohr deformation parameters has been added. 
Restrictions on the complexity of the problem 

The main restriction is the CPU time required for calculations of heavy deformed nuclei 
and for a given precision required. One symmetry plane is assumed. Pairing correlations 
are only included in the BCS limit and for the conserved time-reversal symmetry (i.e., for 
non- rotating states in even-even nuclei). 

Typical running time 

One Hartree-Fock iteration for the super deformed, rotating, parity conserving state of 
1 ggDy 86 takes about nine seconds on the CRAY C-90 computer. Starting from the Woods- 
Saxon wave functions, about fifty iterations are required to obtain the energy converged 
within the precision of about 0.1 keV. In case when every value of the angular velocity 
is converged separately, the complete superdeformed band with precisely determined dy- 
namical moments can be obtained within one hour of CPU on the CRAY C-90, or 
within three hours of CPU on the Athlon-550 MHz processor. This time can be often re- 
duced by a factor of three when a self-consistent solution for a given rotational frequency 
is used as a starting point for a neighboring rotational frequency. 

Unusual features of the program 

The user must have an access to the NAGLIB subroutine f02axe or to the ESSL 
or LAPACK subroutine zhpev which diagonalize complex hermitian matrices, or pro- 
vide another subroutine which can perform such a task. The LAPACK subroutine zh- 
pev can be obtained from the Netlib Repository at University of Tennessee, Knoxville: 
(ittp: / / netlib2.cs.utk.edu/cgi-bin/netlibfiles.pl?filename= /lapack/complex!6 / zhpevl 

LONG WRITE-UP 



1 Introduction 

The method of solving the Hartree-Fock (HF) equations in the Cartesian harmonic oscilla- 
tor (HO) basis was described in the previous publication, Ref. [0, which is below referred 
to as I. The previous version of the code hfodd (vl.60r) was published in Ref. |J which 
is below referred to as II. The present paper is a long write-up of the new version (vl.75r) 
of the code hfodd. This extended version is fully compatible with all previous versions. 

Information provided in I and II remains valid, unless explicitly mentioned in the 
present long write-up. 

In Section |2] we briefly review the modifications introduced in version (vl.75r) of the 
code hfodd. In particular, in Section |2.(j| we present numerical tests of the code in 
situations when the iteration procedure does not converge because of a crossing of single- 
particle levels in function of the rotational frequency. We then introduce a procedure to 
follow configurations along the so-called diabatic path. Such a procedure was implemented 
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in version (vl.75r), and allows, in many cases, to obtain converged solutions for both 
crossing configurations. 

Section |3| lists all additional new input keywords and data values, introduced in version 
(vl.75r). The structure of the input data file remains the same as in the previous versions, 
see Section 3 of II. Similarly, all previously introduced keywords and data values retain 
their validity and meaning. 

2 Modifications introduced in version (vl.75r) 
2.1 Calculation of the density As 

There are two misprints in formulas (1-46)0 for the Cartesian components of the vector 
density As. The correct expressions are as follows: 

As 1 = 25ft + L- + ) +2Ti, (la) 
As 2 = -23 - L-+) + 2T 2 , (lb) 
As 3 = 25? [L ++ - L-) + 2T 3 . (lc) 

The misprints were present only in the text of I, and did not affect the code hfodd 
(vl.60r). Unfortunately, the factors of 2, which should multiply densities T\ and T 2 in 
Eqs. flTa|) and ([TED, respectively, where missing in version (vl.60r) of the code. As far as 
the numerical values are concerned, this error has been fairly unimportant for the final 
results, however, it may have created a weak spurious dependence of the results on the 
orientation of the nucleus with respect to the Cartesian reference frame, because only the 
x and y components of As were affected. 

Incorrect expressions for the densities Asi and As 2 amounted to adding the erroneous 
term — C As (s xt T xt + s yt Ty t ) to the time-odd energy density H° dd (r) of Eq. (I-12a), and 
simultaneously adding the erroneous terms -2C t As T xt and -2C t As T yt to the time-odd spin 
potentials H xt and H yt , respectively, Eq. (I-29b). Therefore, up to a very small difference 
between the matrix elements of spin and kinetic-spin potentials, the incorrect densities 
Asi and As 2 where equivalent to a (direction-dependent) modification of the coupling 
constant Cj. This is the main reason why the error went undetected for a relatively long 
time. Needless to say, calculations performed with C As =0, and in particular those with 
all time-odd terms neglected, are unaffected. 

Since the term s t ■ As t gives anyhow fairly small contribution to the rotational proper- 
ties of nuclei (compare curves denoted by open circles with those denoted by full squares 
in Figs. 3 and 4 of Ref. ||), the incorrect expressions for densities Asi and As 2 had nu- 
merically relatively small importance. The total energies could have been affected at the 
level of about 0.3 MeV and the total spins at the level of about 0.3 h (for details compare 
the output file reproduced in the section TEST RUN OUTPUT below with that given in 

n)- 

3 Symbol (1-46) refers to Eq. (46) of I. 
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2.2 Calculation of the multipole moment Q 2 2 

In version (vl.60r) of the code hfodd, a factor of was missing in the values of the 
multipole moment Q22 printed on the output file. 

Values of the quadrupole moments Q 2 o and Q 22 depend on the normalization factors 
which in the standard nuclear-physics applications are different than those used in elec- 
trodynamics, cf. Ref. In particular, the code hfodd (vl.75r) uses the normalization 
factors such that 

Q 20 = (2z 2 -x 2 -y 2 ), (2a) 
Q22 = Q20 x tan 7, (2b) 

where 7 is the standard angle measuring the non-axial quadrupole deformation ||. 

2.3 Scaling of the coupling constants 

As described in Section 3.2 of II, the code hfodd allows for the scaling of the Skyrme- 
functional coupling constants by arbitrary factors which are read from the input data 
file. The scaling can be independently performed in the total-sum and isoscalar-isovector 
representations of the coupling constants, see Eqs. (1-14) and (1-15). 

In version (vl.60r), the scaling formulas were incorrectly coded and as a result, under 
some very special circumstances described below, calculations were performed for values of 
the scaled coupling constants which were different from those intended by the user. This 
was happening provided that for the given coupling constant two conditions occurred 
simultaneously: 

1. Scaling factor of the isovector coupling constant was different than 1. 

2. Scaling factor of the isovector coupling constant was different than that of the 
isoscalar coupling constant. 

This is a fairly unusual combination of scaling conditions, and hence the error escaped 
detection for a relatively long time. In particular, all results obtained in Ref. || are 
correct because they were obtained with equal scaling factors of the isovector and isoscalar 
coupling constant, and hence condition 2. above was not fulfilled. Obviously, calculations 
performed without scaling, i.e., with all scaling factors equal 1, are unaffected. 

2.4 Diagonalization subroutines 

As described in Section |5.1| , an interface to the LAPACK subroutine zhpev has been 
created. This allows using public-domain diagonalization subroutines, as an alternative 
of using the NAGLIB or ESSL packages, previously described in Section 5.3 of II. 

2.5 Terminating the Hartree-Fock iteration procedure 

As described in Section |3.1| , the code hfodd (vl.75r) can detect several conditions which 
allow to stop the iteration procedure before the predefined number of iterations is ex- 
hausted. This includes stopping after the convergence is reached, and a required stability 
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of the solutions is obtained, as well as an early stopping of the iteration procedure which 
is not going to converge. 

2.6 Diabatic blocking 

Crossing or non-crossing ("repulsion") of the mean-field single-particle energy levels in 
function of some continuous parameter(s) of theory, such as the rotational frequency 
and/or the constraining multipole moments, is one of the most important phenomena in 
the studies of high-spin nuclear states. They may give rise to several kinds of character- 
istic irregularities occurring along the rotational bands and these irregularities often help 
significantly in assigning the theoretical single-particle configurations to the experimental 
rotational bands. 

As it is well known, levels that differ in terms of their discrete good quantum numbers 
such as, for instance, simplex, parity, signature ... etc., will generally cross. These cross- 
ings are such that the underlying intrinsic characteristics (e.g. single-particle alignments 
or multipole moments) do not change in any remarkable way before and after the crossing 
and this mechanism will be of no interest here. 

In contrast, levels that belong to the same discrete symmetry will most of the time 
approach each other in function of the parameter studied and then go apart but in such 
a way that the intrinsic characteristics of the first one will go over to the second one and 
vice versa. This non-crossing rule, sometimes called the Landau-Zener effect, cannot be 
given any more rigorous general formulation, and it may happen that the distance of the 
closest approach for the same symmetry levels is zero. In those, in practice very rare 
cases, the two crossing mechanisms mentioned do not differ. 

An example of a crossing of two HF configurations (no pairing) is presented in Fig. [I]. 
The left panel shows the total energies E(I) as functions of the total spin /, while the 
right panel shows the total Routhians R(u)=E(u)—I(u)u as functions of the rotational 
frequency Tiuj. Since both quantities vary very rapidly as functions of their arguments, 
they are plotted with respect to the corresponding rigid-rotor reference values, i.e., E(I) 
is shown relatively to I(I+1)/(2J ), and R(u) is shown relatively to — J cg> 2 /2, where 
the constant rigid-rotor moment of inertia of Jo=lOOh 2 /MeV has been used. Calcula- 
tions have been performed within the cranking approximation for values of the rotational 
frequency of ^=0.5(0.05)0.8 MeV. 

The examples shown in Fig. [I] correspond to two bands in 151 Tb (see Ref. || for a 
more complete description of calculations performed in this nucleus). Both configurations 
contain the same set of the single-particle levels of 150 Tb being occupied. They correspond 

to the neutron and proton configurations of (N ++ , N + _, N_ + , N ) = (22, 21,21, 21) and 

(15,16,17,17), respectively, as described by the numbers of states occupied in the parity- 
signature blocks (7r,r) = (+l,+i) , (+1,— i), (— 1,+i), and ( The ground-state band 
in 151 Tb can be obtained by putting the 86th neutron into the lowest available level for 
iV + _=22, thus obtaining the closed iV=86 SD magic neutron configuration. Since the 
order of orbitals may change with changing rotational frequency, numbers N n i r are not 
necessarily the most practical for defining physical characteristics of the single-particle 
states in question. Usually, one uses the so-called asymptotic Nilsson quantum numbers 
[Nn z A]Q || for that purpose. Code hfodd calculates these quantum numbers by finding 
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the dominant Nilsson components of the HF single-particle states. Several excited bands 
in 151 Tb can be obtained by putting the 86th neutron into one of the higher available levels, 
for instance, in the iV„ + =22, 23, or 24 levels, that correspond to one of the [521]3/2(r=+i), 
[514]9/2(r=+z), or [761]3/2(r=+i) Nilsson orbitals. 




" I i I i I i L 



50 60 70 80 0.5 0.6 0.7 0.8 
I (fi) -hco (MeV) 

Figure 1: Total energies E(I), (a), and total Routhians R(u), (b), of the [761]3/2(r=-N) 
and [514]9/2(r=+z) neutron diabatic configurations in 151 Tb. The rigid-rotor reference 
energies for Jo=100 h 2 /MeV have been subtracted. (In our graphical representation we 
employ the convention according to which the levels that carry the smoothly varying 
intrinsic characteristics (see text) are denoted with the same symbols. This is at variance 
with the convention that stresses the Landau-Zener mechanism, as used by some authors. 
According to the latter one, the sequences of the lowest-lying points would have been 
drawn as squares and that of the higher- lying points as circles.) 

Of particular interest in 151 Tb are exited bands in which the [761]3/2(r=-N) intruder 
orbital is occupied. This orbital decreases in energy with increasing rotational frequency, 
and therefore, it crosses the other (n,r)= (— orbitals. In particular, in 151 Tb the 
[761]3/2(r=+z) orbital corresponds to the AL + =24th orbital at low frequencies, then it 
becomes the iV_ + =23rd orbital, and finally, at high frequencies it is the lowest available 
iV_ + =22nd orbital. Configurations shown in Fig. [I] correspond to the crossing of the 
[761]3/2(r=+z) and [514]9/2(r=+z) orbitals. Following the standard convention, these 
configurations are called the diabatic ones, because they correspond to the given orbital 
being occupied, irrespective of its excitation energy. On the other hand, configurations 
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based on occupying the AL + =22, 23, or 24 states are called the adiabatic ones. Obviously, 
in adiabatic configurations, different Nilsson orbitals are occupied at different frequencies. 

Fig. |2] shows the negative-parity single-particle neutron Routhians in 151 Tb. The left 
and right panels show the results obtained for the [761]3/2(r=-N) and [514]9/2(r=+i) 
diabatic configurations, respectively. The occupied orbitals are in both cases denoted 
by the filled symbols. It is clear that the crossing frequency depends on which of the 
shown orbitals is occupied. This is due to the self-consistent effects that influence the 
deformations, spins, and other characteristics of many-body states, calculated at given 
rotational frequencies and for given particle-hole configurations. 




Figure 2: Negative-parity neutron single-particle Routhians in 151 Tb calculated for the 
[761]3/2(r=+z) (a) and [514]9/2(r=+z) (b) neutron diabatic configurations. Solid and 
dashed curves denote the r=+i and r=—i signatures, respectively. The arrows denote the 
angular frequencies where the converged solutions near the crossing points could not be 
found. 

In the calculations, one has to begin the analysis by finding the adiabatic configura- 
tions, in order to find the lowest available orbitals, to determine their physical character- 
istics, and to check whether or not there are any crossings of configurations that should be 
followed diabatically. The input data files tbl51-a.dat and tbl51-b.dat give examples 
of runs that find the adiabatic configurations for the iV_ + =23rd and i\L + =24th orbitals 
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being occupied, respectively. 

After several HF iterations, the levels of the same symmetry, that are energetically 
sufficiently far, change their intrinsic characteristics (such as, e.g., angular momentum 
alignments) only by very small amounts. This is not true for the case of the crossing when 
the two states in question exchange their relative positions from one iteration to another. 
Therefore, in the HF calculations, very often the Landau-Zener avoided crossings of levels 
manifest themselves in the form of diverging iteration procedure. This observation is used 
by the code hfodd (vl.75r) in order to locate the crossing by recognizing the oscillatory 
( "ping-pong" ) behavior of solutions. 

Table 1: Example of the output printed when the "ping-pong" divergence is found, see 
text. 

* * 

* TABLE BELOW GIVES THE MAXIMUM VALUES OF DIFFERENCES OF ALIGNMENTS BETWEEN * 

* THE LAST TWO ITERATIONS. IT MAY SERVE AS A GUIDE TO PIN DOWN STATES WHICH * 

* EXCHANGE WAVE FUNCTIONS IN EVERY SECOND ITERATION. THEY MAY BE CANDIDATES * 

* FOR THE DIABATIC BLOCKING CALCULATION * 

* * 

* * * 

* BLOCKS * ANGULAR MOMENTUM ALIGNMENT INTRINSIC SPIN ALIGNMENT * 

* * # 

* * PARTICLES HOLES PARTICLES HOLES * 

* ISO. PAR. SIG. * INDEX VALUE INDEX VALUE INDEX VALUE INDEX VALUE * 

* * * 

* * * 

* NEUT + + * 

* NEUT + - * 

* NEUT * 

* NEUT + * 

* PROT + + * 

* PROT + - * 

* PROT * 

* PROT + * 

* * * 

* * 

* FOLLOWS THE LIST OF CANDIDATE CONFIGURATIONS FOR THE DIABATIC BLOCKING * 

* * 

* * # 

* NEUT - + * 24 2.820 23 2.807 * 

* NEUT + * 24 0.164 23 0.167 * 

* * * 



38 


0, 


.221 


13 


0, 


.018 


24 


0, 


.023 


13 


0, 


.016 


* 


28 





.163 


7 


0, 


,041 


33 


0, 


.014 


7 


0, 


,032 


* 


39 


0. 


.095 


13 


0. 


.198 


44 


0, 


.015 


14 


0. 


.150 


* 


24 


2. 


.820 


23 


2. 


.807 


24 


0, 


.164 


23 


0. 


,167 


* 


24 


0. 


.255 


13 


0, 


.149 


25 


0, 


.129 


14 


0, 


.122 


* 


27 


1. 


.733 


13 





.062 


21 


0, 


.319 


8 


0, 


.020 


* 


32 


0. 


.119 


10 


0. 


.530 


32 


0, 


.029 


10 


0. 


.065 


* 


24 


2 


.649 


15 


0. 


,082 


24 


0, 


.153 


10 


0. 


.030 


* 



The "ping-pong" divergence is a result of the standard self-consistent prescription, 
for which the occupied states in a given iteration are the eigenstates of the mean-field 
Hamiltonian diagonalized in the preceding iteration. Such a divergence is characterized 
by two series of HF states, each one appearing in every second HF iteration. Focusing 
our attention on the specific example discussed above, one can describe such a diver- 
gence in the following way. Suppose that in a given HF iteration, the single-particle 
orbital [761]3/2(r=-M) appears below the [514]9/2(r=-M) orbital, i.e., it corresponds to 
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iV_ + =23. In this particular iteration it becomes occupied, and is included in all densi- 
ties, for which the mean-field Hamiltonian is determined. In the next iteration, the new 
Hamiltonian is diagonalized, however, at present the orbital [514]9/2(r=+i) appears be- 
low [761]3/2(r=+i), and it is this orbital, again for iV_ + =23, which becomes occupied. It 
is obvious that such a situation may repeat itself infinitely many times, never leading to 
a self-consistent solution. 

In many cases, it is possible to find converged solutions near the crossing points simply 
by decreasing the convergence rate, i.e., by using parameters SLOWEV and SLOWOD (see 
Sec. 3.5 of II) only slightly smaller then 1. This method works well provided that the 
two crossing orbitals have fairly different characteristics, like it is the case in the example 
discussed here. However, in the case of a strong mixing between the two crossing orbitals, 
the HF iteration procedure converges near the crossing points only seldomly. 

Upon recognizing the "ping-pong" divergence condition (see Sec. |3.1| for details of the 
procedure), the code hfodd (vl.75r) prints the summary table that may help in identify- 
ing the crossing orbitals. For the iV_ + =23 adiabatic configuration at hu)=0.65 MeV, the 
table has the form presented in Table |I| 

In the upper part of the table, the code prints the absolute values of differences \Ai y \ 
(\As y \) of single-particle average alignments i y (intrinsic spins s y ) between the last and 
the last but one iterations (columns denoted by VALUE) . The results are printed separately 
for the particle (empty) and hole (occupied) states, and their indices are given in columns 
denoted by INDEX. In each of the charge-parity-signature (or charge-simplex) blocks the 
results are printed only for those states for which the values of \Ai y \ and | As y | are the 
largest ones. 

In the lower part of the table, the code prints the list of particle-hole pairs from the 
upper part of the table, for which the indices (INDEX) differ by one, i.e., which correspond 
to the particle-hole pairs at the Fermi surface. In the specific example discussed here, 
the code properly identifies the crossing of the iV_ + =23 and iV_ + =24 orbitals, and pro- 
poses the diabatic blocking based on the average values of the single-particle alignments, 



see Sec. |3.2| . The crossing orbitals, [761]3/2(r=+i) and [514]9/2(r=+z), have markedly 
different alignments of about z y =+2.23 and —0.46 7?,, respectively. This corresponds to 
values of |Aij,| =+2.820 or 2.807 h printed in the table. Therefore, by requesting that the 
state that has a larger alignment be occupied (among the iV_ + =23 and iV_ + =24 orbitals), 
one obtains the diabatic [761]3/2(r=+z) configuration. Similarly, by choosing the state 
that has a smaller alignment, one obtains the [514]9/2(r=+z) diabatic configuration, and 
this is so irrespective of which orbital has smaller or larger energy, i.e., irrespective of the 
angular frequency. The input data files tbl51-c.dat and tbl51-d.dat give examples 
of runs that find the diabatic configurations for the [761]3/2(r=+z) and [514]9/2(r=+z) 
orbitals being occupied, respectively. 

2.7 Saving the Coulomb data 

As described in Section |3.3|, the code hfodd (vl.75r) can save on disk auxiliary data for a 



faster calculation of the Coulomb potential. These data can be reused in subsequent runs 
of the code, provided the Coulomb parameters NUMCOU, NUMETA, and FURMAX (see Section 
3.5 of II) do not change from one run to the next. 
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2.8 Calculation of moments and radii 



Average quadrupole moments, (920) and (922), radii squared, (r 2 ), and sizes squared, (x 2 ), 
(y 2 ), and (z 2 ), of single-particle states are calculated and stored in the review file FILREV, 
see Section 3.9 of II. 

2.9 Calculation of the Bohr deformation parameters 

In version (vl.75r), the Bohr deformation parameters q.\^ |5| are calculated using the first- 
order approximation. The code uses the simplest linearized expression (cf. Eq. (1.35) of 
Ref. ||) relating deformation parameters to multipole moments of a sharp-edge uniformly 
charged shape, i.e., 



For /i^O the deformation parameters contain the standard factor of This reflects the 
fact that for the assumed S y simplex symmetry (see Sec. 3.1 of I), the values of multipole 
moments for /x<0 are up to a phase equal to those for /x>0. Neutron, proton, or mass 
deformation parameters are calculated from the corresponding neutron, proton, or mass 
multipole moments (r A Y\ M ), and by setting N in Eq. (Q) equal to the number of neutrons, 
protons, or nucleons. The equivalent radii Rq are respectively calculated from the neutron, 
proton, or mass rms radii as ^/5/3(r 2 ) 1//2 . One should note that for large deformations the 
neglected higher-order terms (see Eq. (1.35) of Ref. ||) will in general be non- negligible. 
Therefore, the printed values of the Bohr deformation parameters should be used only as 
(often very) rough estimates. 

3 Input data file 

Structure of the input data file has been described in II, and in version (vl.75r) of the 
code hfodd this structure is exactly the same. All previous items of the input data file 
remain valid, and several new items were added, as described in Sees. |3.1| - 0.3| . 



Together with the FORTRAN source code in the file hfodd. f, several examples of the 
input data files are provided. File dyl52-c.dat contains all the valid input items, and 
the input data are identical to the default values. Therefore, the results of running the 
code with the input data file dyl52-c.dat are identical to those obtained for the input 
data file containing only one line with the keyword EXECUTE. 

File dyl52-d.dat contains all those input items that were added in version (vl.75r), 
and the input data there are those recommended for being used in the new version. This 
file is reproduced in the section TEST RUN INPUT below. 

3.1 General data 
Keyword: ITERAT_EPS 

0.0 = EPSITE 

Keyword: MAXANTIOSC 

1 = NULAST 



11 



These two parameters govern the termination of the HF iterations according to the 
achieved stability of solutions. The stability of the HF energy has been defined in (1-37) 
as the difference between the total energies calculated from the single-particle energies 
and from the Skyrme functional. The HF iterations continue until the absolute value of 
the stability is smaller than EPSITE (in MeV) over NULAST consecutive iterations. When 
this condition is fulfilled, iteration procedure terminates and the final results are printed. 
This allows for an automated adjustment of the number of iterations which are required 
to achieve a given level of convergence. The number of iterations NOITER, Section 3.1 of 
II, can now be set to a large value at which the iterations terminate if a stable solution is 
not found. 

The default value of EPSITE=0.0 ensures that whenever this new option is not used, 
the code hfodd (vl.75r) behaves as that in version (vl.60r). If a non-zero value of EPSITE 
is used, a non-zero value of NULAST has to be used too. In practice, a value of NULAST=5 
prevents an accidental termination of iterations in all cases when the stability energy 
changes the sign, but the solution is not yet self-consistent. 

Keyword: PING-PONG 

0.0, 3 = EPSPNG, NUPING 
The code is able to detect the "ping-pong" divergence described in Section |2]6|, i.e., 
the situation when the HF iteration procedure gives alternating results in every second 
iteration. Upon continuing the iteration, both sequences of results, i.e., those which 
correspond to the iteration numbers being even and odd, stay different but perfectly 
stable, and hence the correct self-consistent solution is never attained. 

The code recognizes such a situation by calculating the averages and variances of the 
stability energy (1-37), separately in the even and in the odd sequences of results, over 
the last NUPING pairs of iterations. The "ping-pong" divergence condition occurs when 
both variances become a factor EPSPNG smaller than the absolute value of the difference 
of the corresponding averages, i.e., when 



where 



and 



A(5£) cvcn < EPSPNG x \5£ even - 5£ odd \ 
A(SS) odd < EPSPNG x \6S even — ~5£ odd \ 



(NUPING-1 
E 5£ n-2p ] /NUPING, 

(NUPING-1 
]T 5£ n - 2p -i ) /NUPING, 
p=0 



A(5£)even 

A(5S) odd 



(NUPING-1 \ 1//2 

E ( S£ n-2p ~ SEeven) J /NUPING, 

/nuping-i . . 2 ^ 1//2 



I E (^n-2p-l _ ^£odd 
\ P=0 



/NUPING. 



(4a) 
(4b) 



(5a) 
(5b) 

(6a) 
(6b) 
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Here, n denotes the number of the last accomplished HF iteration, and 5S n / denotes the 
stability energy (1-37) obtained in the n'-th iteration. 

The default value of EPSPNG=0.0 ensures that whenever this new option is not used, 
the code hfodd (vl.75r) behaves as that in version (vl.60r). If a non-zero value of EPSPNG 
is used, a value of NUPINOl has to be used too. In practice, values of EPSPNG=0.01 and 
NUPING=3 allow for an efficient detection of the "ping-pong" divergence condition. 

Upon discovering the "ping-pong" divergence, the HF iterations are terminated and 
a table of absolute values of maximum differences of single-particle observables between 
the two sequences of iterations is printed, see Sec. and Table [TJ. These maximum 
differences are determined for states in each of the charge-parity-signature (or charge- 
simplex) blocks, and separately for particle and hole states. Whenever such a maximum 
difference is found for a particle state and for a hole state with adjacent indices, such a 



pair is proposed as a candidate for the diabatic blocking calculation, see Sees. |2.6| and |372 
Keyword: CHAOTIC 

= NUCHAO 

The code is able to detect the "chaotic" divergence which occurs when the HF iterations 
give results which chaotically vary from one iteration to another one. The code recognizes 
such a divergence by finding the local maxima M^, k—1,2,. . . , in the sequence of absolute 
values of the stability energies (1-37), obtained in the entire series of the HF iterations 
performed. The "chaotic" divergence condition occurs when the code finds NUCHAO positive 
differences Mj.-Mj._i. When this condition occurs, iteration procedure terminates and 
the final results are printed. 

For NUCHA0=0 (the default value) the code does not check whether the "chaotic" 
divergence occurs or not. In practice, a value of NUCHA0=5 allows for an efficient detection 
of the "chaotic" divergence condition. However, for a small value of NUCHAO and a small 
value of EPSPNG, the "ping-pong" divergence can sometimes be mistaken for the "chaotic" 
divergence. If one is interested in the diabatic-blocking data, printed after the "ping-pong" 
divergence, the recommended value of NUCHA0=5 should be increased to 10 or more. 

Keyword: PHASESPACE 

0, 0, 0, =NUMBSP(0,0), NUMBSP (1,0), 
NUMBSP(0,1), NUMBSP (1,1) 

Numbers of the lowest mean-field eigenstates which are kept after the diagonalization of 
the mean-field Hamiltonians in the four simplex-charge blocks: (s,q) = (+i,n), (— i, n), 
(+i,p), (—i,p). All other eigenstates are discarded. If any of these numbers is equal to 
zero (the default value), the code sets it equal to the number of neutrons IN_FIX (for q=n) 
or protons IZ_FIX (for q=p), see Section 3.1 of II. 

For calculations without pairing, the user is responsible for using values of NUMBSP 
large enough to accommodate all wave functions which might be useful for the required 
vacuum and particle-hole configurations, see Section 3.4 of II. In practice, the use of the 
default values described above is recommended as a safe option. The execution time is 
almost independent of NUMBSP. The size of the matrices defined by the NDSTAT parameter 
can be reduced by the user for smaller values of NUMBSP. 
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3.2 Configurations 



Keyword: DIASIM_NEU 

2,2, 1,1, 0,0= KPFLIM(0,0),KPFLIM(1,0), 

KHFLIM (0 , 0) ,KHFLIM (1,0), 
KOFLIM (0 , 0) ,K0FLIM (1,0) 
Diabatic blocking of neutron single-particle simplex configurations. Matrices KPFLIM con- 
tain the indices of the particle states in the two simplex blocks denoted by (+) and (— ), of 
given simplex values, i.e., s—+i and —i, respectively. Matrices KHFLIM contain analogous 
indices of the hole states, and matrices KOFLIM define type of blocking according to the 
following table: 



K0FLIM=0 •<=>- No diabatic blocking in the given parity/signature block. 
K0FLIM=+1 •<=>- The state which has the larger alignment is occupied. 
K0FLIM=— 1 •<=>- The state which has the smaller alignment is occupied. 
K0FLIM=+2 •<=>- The state which has the larger intrinsic spin is occupied. 
K0FLIM=— 2 •<=>- The state which has the smaller intrinsic spin is occupied. 



Within the diabatic blocking procedure one does not predefine whether the particle or 
the hole state is occupied (like is the case when the particle-hole excitations are defined, see 
Section 3.4 of II). In each iteration the code calculates the average alignments (or average 
intrinsic spins) of both states (those defined by KPFLIM and KHFLIM), and occupies that 
state for which a larger, or a smaller value is obtained. Therefore, the order of both states 
in the Routhian spectrum is irrelevant. 

The user is responsible for choosing the particle-state indices (in KPFLIM) only among 
those corresponding to empty single-particle states, and the hole-state indices (in KHFLIM) 
only among those corresponding to occupied single-particle states, see Section 3.4 of II. 

Keyword: DIASIM_PR0 

2,2, 1,1, 0,0= KPFLIM(0,1),KPFLIM(1,1), 

KHFLIM (0,1) ,KHFLIM (1,1), 
K0FLIM(0,1),K0FLIM(1,1) 

Same as above but for the diabatic blocking of proton single-particle simplex configura- 
tions. 

Keyword: DIASIG_NEU 

2,2,2,2, 1,1,1,1, 0,0,0,0 = 

KPFLIG (0 ,0,0) ,KPFLIG (0 ,1,0) ,KPFLIG (1,0,0) ,KPFLIG (1,1,0), 
KHFLIG (0 ,0,0) ,KHFLIG (0 ,1,0) ,KHFLIG (1,0,0) ,KHFLIG (1,1,0), 
K0FLIG(0,0,0),K0FLIG(0,1,0),K0FLIG(1,0,0),K0FLIG(1,1,0) 
Diabatic blocking of neutron single-particle parity/signature configurations. Matrices 
KPFLIG contain the indices of particle states in the four parity/signature blocks de- 
noted by (+,+), (+,—), (— ,+), and (—,—), of given (paritysignature) combinations, i.e., 
(7r, r)—(+l,+i), and respectively. Matrices KHFLIG contain 

analogous indices of hole states, and matrices K0FLIG define the type of blocking accord- 
ing to the table of values identical to that defined for the simplex case above. Other rules 
described for the simplex case apply here analogously. 
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Keyword: DIASIG_PRO 

2, 2, 2, 2, 1, 1, 1, 1, 0, 0, 0, = 

KPFLIG (0 ,0,1) ,KPFLIG (0 ,1,1) ,KPFLIG (1,0,1) ,KPFLIG (1,1,1), 
KHFLIG (0 ,0,1) ,KHFLIG (0 ,1,1) ,KHFLIG (1,0,1) ,KHFLIG (1,1,1), 
K0FLIG(0,0,1),K0FLIG(0,1,1),K0FLIG(1,0,1),K0FLIG(1,1,1) 

Same as above but for the diabatic blocking of proton single-particle parity/signature 
configurations. 

3.3 Files 
Keyword: C0UL0MFILE 

HFODD.COU = FILC0U 
CHARACTER*68 file name of the file containing auxiliary data for a faster calculation of the 
Coulomb potential. Must start at the 13-th column of the data line. Parameters NUMC0U, 
NUMETA, and FURMAX which define the calculations of the direct Coulomb potential (see 
Section 5 of I and Section 3.5 of II) are usually kept unchanged for the whole series of 
calculations in one region of nuclei. Therefore, many Coulomb auxiliary results can be 
calculated only once, and stored in the file FILC0U. Handling of this file is determined by 
the input parameters IC0ULI and IC0UL0. 

Keyword: C0UL0MSAVE 

0,0 = IC0ULI, IC0UL0 

Input parameters IC0ULI and IC0UL0 determine actions pertaining to reading and/or 
writing of the Coulomb file FILC0U, according to the following table: 

IC0ULI IC0UL0 Action 

neither read nor write the Coulomb file 

1 read, but do not write the Coulomb file 

1 do not read, but write the Coulomb file 

1 1 use automated handling of the Coulomb file 

The default values of IC0ULI=IC0UL0=0 ensure that whenever this new option is not 
used, the code hfodd (vl.75r) behaves as that in version (vl.60r). However, unless it 
is required by special circumstances, a use of the automated handling of the Coulomb 
file (IC0ULI=IC0UL0=1) is recommended. Within the automated mode, the code checks 
whether the Coulomb file exists, and whether it contains data which match the current 
values of the input parameters NUMC0U, NUMETA, and FURMAX. If this is the case, the code 
reads the data from the Coulomb file. If this is not the case, the code calculates the 
Coulomb auxiliary results and stores them in the Coulomb file FILC0U. In the automated 
mode, the user is informed by appropriate messages printed on the output file about what 
type of the action has been taken in the given run of the code. 

4 Output file 

Together with the FORTRAN source code in the file hfodd. f, an example of the output 
file is provided in dyl52-d.out. Selected lines from this file are presented in the section 
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TEST RUN OUTPUT below. This output file corresponds to the input file dyl52-d.dat 
reproduced in the section TEST RUN INPUT below. 



5 FORTRAN source file 

The FORTRAN source code in is provided in the file hfodd.f and can be modified in 
several places, as described in this section. 

5.1 Library subroutines 

The code hfodd requires an external subroutine which diagonalizes complex hermitian 
matrices. Version (vl.75r) contains an interface to the LAPACK subroutine zhpev that 
can be downloaded (with dependencies) from 

jittp: / /netlib2.cs.utk.edu/cgi-bin/netlibfiles.pr.''filename=/lapack/complexl6/zhpevl 
This subroutine and its dependencies are in the REAL*8 and COMPLEX* 16 versions, and 
should be compiled without promoting real numbers to the double precision. On the other 
hand, the code hfodd itself does require compilation with an option promoting to double 
precision. Therefore, the code and the zhpev package should be compiled separately, and 
then should be linked together. 

In order to activate the interface to the LAPACK zhpev subroutine, the following 
modifications of the code hfodd (vl.75r) have to be made: 

1. Change everywhere the value of parameter I_CRAY=1 into I_CRAY=0. 

2. Change everywhere the value of parameter IZHPEV=0 into IZHPEV=1. 

3. Change the name of the subroutine zhpev, provided with the code hfodd, to 
another name, or remove it from the file. 

4. If your compiler does not support undefined externals, or subroutines called with 
different parameters, remove calls to subroutines cgemm, f02axe, and the first call 

to ZHPEV. 
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TEST RUN INPUT 



I This file (dyl52-d.dat) contains the input data for the code HFODD. I 
I Only keywords introduced after version (1.60r) are included here. I 



General data 

ITERAT_EPS 

0.0001 

MAXANTIOSC 

5 

PING-PONG 

0.01 3 

CHAOTIC 



PHASESPACE 



Configurations 



DIASIM. 


_NEU 




PARTICLES 
00 00 






HOLES 
00 00 






TYPE 





DIASIM. 


.PRO 




PARTICLES 
00 00 






HOLES 
00 00 






TYPE 





DIASIG. 


_NEU 




PARTICLES 






HOLES 






TYPE 








00 


00 00 


00 


00 


00 00 


00 











DIASIG. 


.PRO 




PARTICLES 






HOLES 






TYPE 








00 


00 00 


00 


00 

- Files 


00 00 


00 












COULOMFILE 

HFODD . COU 

COULOMSAVE 

1 1 

Calculate 

EXECUTE 
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TEST RUN OUTPUT 



******************************************************************* 

* HFODD HFODD HFODD HFODD HFODD HFODD HFODD HFODD * 

* * 
******************************************************************************* 

* * 

* SKYRME-HARTREE-FOCK CODE VERSION 1.75R * 

* * 

* ONE SYMMETRY-PLANE AND NO TIME-REVERSAL SYMMETRY * 

* DEFORMED CARTESIAN HARMONIC-OSCILLATOR BASIS * 

* * 
******************************************************************************* 

* * 

* JACEK DOBACZEWSKI AND JERZY DUDEK * 

* * 

* INSTITUT DE RECHERCHES SUBATOMIQUES * 

* * 

* (FORMER CENTRE DE RECHERCHES NUCLEAIRES) * 

* STRASBOURG, 1993-2000 * 

* .- .. * 

******************************************************************************* 

* * 

* PARAMETER SET SKM* : T0= -2645.00 Tl= 410.00 T2= -135.00 T3= 15595.00 * 

* P0WER=0.1667 W=130 X0= 0.09000 Xl= 0.00000 X2= 0.00000 X3= 0.00000 * 

* * 
******************************************************************************* 

******************************************************************************* 

* AUTOMATED HANDLING OF THE COULOMB FILE * 

* * 
******************************************************************************* 

* * 

* THE COULOMB FILE WRITTEN ON THE DISC, FILE NAME (FILCOU) : * 

* HFODD . COU * 

* * 
******************************************************************************* 

* NUCLIDE: N = 86 Z = 66 NUMBER OF ITERATIONS = 50 CONTINUATION? = * 

* * 
******************************************************************************* 

* ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITEO . 000100 * 

* FOR THE 5 CONSECUTIVE ITERATIONS * 

* * 

* * 

* ITERATIONS WILL STOP WHEN THE PING-PONG DIVERGENCE OCCURS FOR EPSPNG=0.010 * 

* DURING 3 PAIRS OF CONSECUTIVE ITERATIONS * 

* * 
******************************************************************************* 

* ITERATIONS WILL STOP WHEN THE CHAOTIC-DIVERGENCE CONDITION OCCURS 5 TIMES * 

* * 
******************************************************************************* 

* * 

* SLOW-DOWN PARAMETERS: TIME-EVEN FIELDS =0.50 TIME-ODD FIELDS =0.50 * 

* * 
******************************************************************************* 

* COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * 

* * 

t MAXIMUM NUMBER OF MULTIPOLES CONSIDERED: FOR THE CONSTRAINTS, NMUCON = 2 * 

* FOR THE COULOMB FIELD, NMUCOU = 4 * 

* FOR THE OUTPUT INFO, NMUPRI = 4 * 

* * 

* * 

* MULTIPOLE CONSTRAINTS: LAMBDA= 2 MIU= MOMENT= 42.000 STIFFNESS=0 . 010 * 

* * 
******************************************************************************* 

* * 

* PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST (1) * 

* THE MIDDLE (0) * 

* AND/OR THE LAST (1) * 

* * 
******************************************************************************* 



* NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS 

* NEUTRONS 

* PROTONS 



TIME-UP TIME-DOWN 
86 86 
66 66 



^*)|t^>|E:4;^*;+i:+:^^>|E:4:^:+:>|E:+:)|E:+:>|E:4;^)^ 

* CALCULATIONS WITHOUT PAIRING * 

* * 

* * 

* CALCULATIONS WITH PARITY/SIGNATURE SYMMETRY * 

* * 
******************************************************************************* 

* CALCULATIONS WITH BROKEN TIME-REVERSAL SYMMETRY * 

* * 
******************************************************************************* 

* INITIAL ROTATIONAL FREQUENCY OMEGAY = 0.500000 * 

* * 
******************************************************************************* 

* LINEAR CONSTRAINT ON SPIN * 

* , , ,- . ,- ,- ,- * 
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******************************************************************** 

* * 

* PARITY/SIGNATURE CONFIGURATIONS : * 
* 

* V_A_C_U_U_M P_0_T_I_C_L_E_S 

* (++)~(+ : r( : +)~(— ) (++r(+ : )~( : +r(--) 



H_0_L_E_S 

(++) ~(+ z r~( : +) (~) 



NEUTRONS : 
PROTONS : 



22 
16 



22 
16 



21 
17 



21 
17 



******************************************************************************* 
******************************************************************************* 

* CONVERGENCE REPORT * 

* * 
******************************************************************************* 



* 


ITER 


ENERGY 


STABILITY 


020 


GAMMA 


SPIN 


OMEGA 


HOW NICE 


* 


* 
* 





-obo 


rioo/t f\~7 
. U2y4U / 


-/bl . 


. booo f 1 


54 . 


. 1 / o 


-u . 


069 


yo . 


.III 


u . 


, 500 


u 


. 4zzo4z 


* 

* 


* 


1 


-1163 . 


. o ( zloy 


87 . 


. t 2b»bo 


48 


. 776 


o . 


. 047 


52 . 


, 033 


. 


, 500 


l 


. 081557 


* 


* 


2 


-1195 . 


. yoouoo 


63 . 




47 


. 722 


-0 . 


. 025 


49 . 


, 155 


o . 


, 500 


l 


. ubbyy^ 


* 


* 


3 


-1204 




38 


471958 


46 


175 


o 


003 


45 


243 


o 


500 


i 




* 


* 


4 


-1207! 


. 038466 


24 : 


1979693 


45. 


:544 


o'. 


.020 


45' 


!509 


o' 


!500 


l 


,021132 


* 


* 


5 


-1208. 


.304258 


15 


.158481 


44 


.918 


0. 


.038 


45. 


,757 


0. 


,500 


l 


.012705 


* 


* 


6 


-1208. 


.924216 


9 


.559571 


44 


.422 


0. 


.055 


46. 


,049 


0. 


,500 


l 


.007971 


* 


* 


7 


-1209. 


. 227387 


6. 


.030570 


44. 


.007 


0. 


071 


46. 


,339 


0. 


,500 


l 


,005012 


* 


* 


8 


-1209. 


.364247 


3. 


.884991 


43 


.664 


0. 


086 


46. 


,612 


0. 


,500 


l 


,003223 


* 


* 


9 


-1209. 


.413569 


2. 


.543565 


43 


.379 


0. 


101 


46. 


,858 


0. 


,500 


l 


,002108 


* 


* 


10 


-1209. 


.417195 


1. 


.696859 


43. 


.142 


0. 


.114 


47. 


,074 


0. 


,500 


l 


.001405 


* 


* 


25 


-1209. 


.092516 


0. 


.016427 


42. 


.016 


0. 


.215 


48. 


,199 


0. 


,500 


l 


.000014 


* 


* 


26 


-1209. 


.085246 


0. 


.012900 


42. 


.002 


0. 


.217 


48. 


,214 


0. 


,500 


l 


.000011 


* 


* 


27 


-1209. 


.079025 


0. 


.010203 


41. 


,990 


0. 


220 


48. 


,227 





,500 


l 


,000008 


* 


* 


28 


-1209. 


.073706 


0. 


.008125 


41 


.980 


0. 


.222 


48. 


,238 


0. 


,500 


l. 


.000007 


* 


* 


29 


-1209. 


.069161 





.006513 


41 


.972 


0. 


.223 


48. 


,248 


0. 


,500 


l 


.000005 


* 


* 


30 


-1209 


.065279 





.005253 


41. 


,965 


0. 


225 


48. 


,256 


0. 


,500 


l 


,000004 


* 


* 


31 


-1209. 


.061965 





.004261 


41. 


,959 


0. 


226 


48. 


,263 


0. 


,500 


l 


,000004 


* 


* 


32 


-1209 


.059136 





.003474 


41. 


,954 


0. 


.228 


48. 


,268 


0. 


,500 


l 


,000003 


* 


* 


33 


-1209. 


.056722 





. 002846 


41 


.950 


0. 


.229 


48. 


,273 


0. 


,500 


l 


.000002 


* 


* 


34 


-1209. 


.054662 





.002342 


41 


.946 


0. 


.230 


48. 


,278 


0. 


,500 


l 


.000002 


* 


* 


35 


-1209. 


.052904 





.001934 


41 


.943 


0. 


.231 


48. 


,281 


0. 


,500 


l 


.000002 


* 


* 


36 


-1209. 


.051404 





.001603 


41. 


,941 


0. 


231 


48. 


,284 


0. 


,500 


l 


,000001 


* 


* 


37 


-1209. 


.050124 





.001333 


41 


.938 


0. 


.232 


48. 


,287 


0. 


,500 


l 


.000001 


* 


* 


38 


-1209 


.049031 


0. 


.001112 


41. 


.937 


0. 


.233 


48. 


,289 


0. 


,500 


l 


.000001 


* 


* 


39 


-1209. 


. 048098 





.000929 


41. 


.935 


0. 


.233 


48. 


,291 


0. 


,500 


l 


,000001 


* 


* 


40 


-1209. 


.047301 





.000778 


41. 


.934 


0. 


.234 


48. 


,292 


0. 


,500 


l 


,000001 


* 


* 


41 


-1209. 


.046621 





.000653 


41. 


,933 


0. 


234 


48. 


,294 


0. 


,500 


l 


,000001 


* 


* 


42 


-1209. 


. 046040 


0. 


.000549 


41. 


.932 


0. 


.235 


48. 


,295 


0. 


,500 


l 


.000000 


* 


* 


43 


-1209. 


. 045544 





.000462 


41. 


.931 


0. 


.235 


48. 


,296 


0. 


,500 


l 


.000000 


* 


* 


44 


-1209. 


.045120 





.000389 


41. 


.931 


0. 


.235 


48. 


,297 


0. 


,500 


l 


.000000 


* 


* 


45 


-1209. 


.044757 


0. 


.000329 


41. 


,930 


0. 


235 


48. 


,298 





,500 


l 


,000000 


* 


* 


46 


-1209. 


. 044447 





.000278 


41. 


.930 


0. 


.236 


48. 


,298 


0. 


,500 


l 


.000000 


* 


* 


47 


-1209. 


.044182 





.000235 


41. 


,929 


0. 


.236 


48. 


,299 


0. 


,500 


l 


.000000 


* 


* 


48 


-1209 


043955 





.000199 


41. 


.929 


0. 


236 


48. 


.299 


0. 


,500 


l 


,000000 


* 


* 


49 


-1209. 


.043761 





.000168 


41. 


,929 


0. 


236 


48. 


,300 


0. 


,500 


l 


,000000 


* 



****************************************** 

* DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * 

* * 



* 






TOTAL (T) 




SUM(S) 


ISOSCALAR(P) 


ISOVECTOR(M) 


* 
* 


DRH0 




17.558301 


8. 


. 908547 


17.558301 


0.258793 


* 


DRH0D 




12.481232 


6. 


331119 


12.481232 


0.181006 


* 


DLPR 




-3.718467 


-1. 


.885747 


-3.718467 


-0.053027 


* 


DTAU 




15.576998 


7. 


972052 


15.576998 


0.367107 


* 


DSCU 




0.119021 


0. 


.063593 


0.119021 


0.008164 


* 


DDIV_ 




0.821038 


0. 


.417814 


0.821038 


0.014590 


* 
* 


DSP I 




. 009448 


0. 


005856 


. 009448 


0.002265 


* 


DSPID 




0.006720 


0. 


.004158 


0.006720 


0.001595 


* 


DLPS 




-0.004751 


-0. 


006086 


-0.004751 


-0.007421 


* 


DCUR 




0.058252 


0. 


030270 


0.058252 


0.002287 


* 


DKIS 




0.014893 


0. 


. 007474 


0.014893 


0.000055 


* 


DR0T 




0.004337 


0. 


.002462 


. 004337 


0.000587 



******************************************************************************* 

******************************************************************************* 

* * 

* CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * 

* * 
******************************************************************************* 



* 




TOTAL (T) 


SUM(S) 


ISOSCALAR(P) 


ISOVECTOR(M) 


* 
* 


ERHO = 


-24265. 


791099 


6951. 


.116415 


-17415. 


.639545 


100, 


.964861 


* 


ERHOD = 


16220. 


.401153 


-4113. 


908512 


12165, 


.300865 


-58, 


.808223 


* 


ELPR = 


317. 


.231713 


-64. 


.527908 


253, 


.611067 


-0. 


.907262 


* 


ETAU = 


1070. 


.918626 


-543. 


096073 


540, 


.327125 


-12, 


.504572 


* 


ESCU = 


0. 


.000000 


4. 


.332261 


4. 


.054159 


0. 


.278102 


* 


EDIV_ = 


-53. 


.367500 


-27. 


157931 


-80, 


.051250 


-0. 


.474181 


* 
* 
* 


SUM EVEN: 


-6710. 


.607107 


2206. 


.758253 


-4532 , 


.397579 


28. 


.548725 


* 


ESPI = 


-0. 


.562245 


3. 


.872319 


2. 


.561340 





,748734 


* 


ESPID = 


0. 


.000000 


-2. 


701576 


-2 


. 183390 


-0. 


,518186 


* 


ELPS = 


0. 


.000000 


-0. 


.208253 


-0 


.081291 


-0 


,126962 


* 


ECUR = 


-4. 


.004811 


2. 


062118 


-2 


.020609 


0. 


,077916 


* 


EKIS = 


0. 


.000000 


-0. 


.509167 


-0. 


.507277 


-0. 


,001890 


* 


ER0T_ = 


-0. 


.281919 


-0. 


. 160046 


-0 


.422878 


-0 


.019087 


* 
* 


SUM ODD: 


-4. 


. 848975 


2. 


.355395 


-2. 


.654105 


0. 


.160525 



******************************************************************************* 
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9(C 9|C =|C 9|C =|C 9|C 3|C 3|C 3|C 3|C 3|C =#= >|C 3|C 3|C =#= =t= 3|C =t= =t= ^fc =f= 9(1: 3(C 3f= ^ H= 3^ H= 3)d 3|C 9(C 3^ 9(C 3f= ^ 3fc 9)= 3fc 9)= =fc 9(1: =f= afc =f= afc 3|C =f= 3|C Sfc 

* SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK 
* 

******************************************** 



***************** 

NEUTRONS * 
* 

***************** 



* 


NO) 


ENERGY 


(++ 


+- 


-+ 


— ) 


N 


nz 


lz 


OMEG> 


<P> 




JY 




SY 


GFACT 




7R1 


-11 


812 


( 





20 


01 


7 







l/2> 


-100 


3 


125 





328 





105 


* 


7fi1 


-11 


744 


( 








20) 


5 


3 


2 


3/2> 


-100 





297 


-0 


063 


-0 


211 


* 


771 


-11 


742 


( 





21 


01 


5 


3 


2 


3/2> 


-100 




567 





110 





071 


* 


78") 
f 0/ 


-11 


527 


(19 








01 


4 


I 




1 /9> 
j./ z^ 


100 





130 


-0 


165 


-1 


268 


* 


791 


-11 


320 


( 


19 





01 


4 


I 


3 


5/2> 


100 


-0 


152 


-0 


004 





026 




80) 


-11 


284 


(20 








0) 


4 


1 


3 


5/2> 


100 





169 


-0 


109 


-0 


648 


* 


81 1 


-11 


251 


(21 








01 


g 


5 




l/2> 


100 





771 


-0 


167 


-0 


217 


* 


82) 


-11 


034 


( 


20 





01 


4 


I 




l/2> 


100 





522 


-0 


004 


-0 


008 


* 


83) 


-10 


899 


( 


21 





01 


4 






l/2> 


100 





345 


-0 


016 


-0 


046 


* 


R41 


-10 


416 


(22 








01 


g 


4 


2 


5/2> 


100 


-0 


014 


-0 


031 


2 


177 


* 


85) 


-10 


405 


( 


22 





01 


g 


4 


2 


0/ z^ 


100 





074 


-0 


019 


-0 


263 


* 


86) 


-9 


561 


( 








21 1 


7 


7 





l/2> 


-100 


2 


188 


-0 


036 


-0 


016 


* 


871 


-7 


830 


( 





22 


01 


5 


2 




3/2> 


-100 





839 





124 





148 


* 


88) 


-7 


741 


( 








22) 


5 


2 




3/2> 


-100 





308 





127 





412 


* 


89) 


-7 


677 


(23 








01 


4 





2 


5/2> 


100 


-0 


282 





087 


-0 


308 


* 


90) 


-7 


675 


( 


23 





0) 


4 





2 


5/2> 


100 


-0 


282 





086 


-0 


304 


* 




-7 


521 


( 





23 


01 


5 


2 




3/2> 


-100 





975 





087 





089 


* 


92) 


-7 


215 


( 








23) 


7 


7 





l/2> 


-100 





026 


-0 


124 


-4 


839 


* 


93) 


-7 


208 


( 








941 


5 




4 


9/2> 


-100 


-0 


334 





008 


-0 


024 


* 


941 


-7 


208 


( 





24 


01 


5 




4 


9/2> 


-100 


-0 


336 





009 


-0 


026 




95) 


-6 


936 


( 


24 





01 


g 


4 





l/2> 


100 




271 





325 





256 




96) 


-6 


423 


(24 








01 


g 


4 





l/2> 


100 





341 


-0 


415 


-1 


220 




971 


-6 


364 


( 








25) 


5 


2 


I 


1 /9> 
j./ z^ 


-100 





258 


-0 


215 


-0 


832 




98) 


-6 


263 


(25 








01 


g 


3 


3 


7/2> 


100 


-0 


015 





009 


-0 


602 


* 


99) 


-6 


263 


( 


25 





0) 


g 


3 


3 


7/2> 


100 


-0 


009 





019 


-2 


112 




100) 


-6 


164 


( 


26 





0) 


4 








l/2> 


100 





141 





448 


3 


168 




101) 


-5 


819 


( 





25 


01 


5 


2 




l/2> 


-100 


-0 


045 


-0 


064 


1 


442 




102) 


-5 


703 


(26 








01 


4 








1 /9> 
±/ z^ 


100 


-0 


575 


-0 


399 





694 




103) 


-5 


505 


( 





26 


01 


5 


2 


3 


5/2> 


-100 





095 


-0 


021 


-0 


223 




104) 


-5 


503 


( 








26) 


5 


2 


3 


5/2> 


-100 





074 


-0 


012 


-0 


161 


* 


105) 


-5 


369 


( 


27 





0) 


g 


4 


2 


3/2> 


100 





350 





004 





Oil 


* 


106) 


-5 


366 


(27 








0) 


6 


4 


2 


3/2> 


100 





608 


-0 


019 


-0 


032 


* 


107) 


-5 


172 


( 


28 





0) 


4 





2 


3/2> 


100 


-0 


519 


-0 


021 





041 


* 


108) 


-5 


089 


(28 








0) 


4 





2 


3/2> 


100 


-0 


445 


-0 


152 





341 


* 


109) 


-4 


928 


( 


29 





0) 


4 





4 


7/2> 


100 


-0 


387 


-0 


031 





080 


* 


110) 


-4 


926 


(29 








0) 


4 





4 


7/2> 


100 


-0 


370 


-0 


040 





109 


* 


111) 


-4 


793 


( 


30 





0) 


8 


8 





l/2> 


100 


5 


270 





469 





089 


* 


112) 


-4 


385 


( 








27) 


7 


6 


1 


l/2> 


-100 


1 


575 


-0 


105 


-0 


067 


* 


113) 


-4 


348 


( 





27 


0) 


7 


7 





l/2> 


-100 


1 


733 





235 





136 


* 


114) 


-3 


967 


( 





28 


0) 


7 


5 


2 


5/2> 


-100 





234 


-0 


038 


-0 


161 


* 


115) 


-3 


927 


( 








28) 


7 


5 


2 


5/2> 


-100 


-0 


157 


-0 


155 





985 


* 


116) 


-2 


555 


(30 








0) 


8 


7 


1 


3/2> 


100 


3 


047 





077 





025 


* 


117) 


-2 


480 


( 








29) 


5 


1, 


2 


5/2> 


-100 


-0 


074 





067 


-0 


915 


* 


118) 


-2 


478 


( 





29 


0) 


5 


1, 


2 


5/2> 


-100 


-0 


067 





069 


-1 


033 


* 


119) 


-2 


222 


( 








30) 


5 





5 


ll/2> 


-100 


-0 


602 


-0 


002 





003 



****************************************************************** 



* MULTIPOLE MOMENTS IN UNITS OF (10 FERMI) **LAMBDA 
* 



* 

TOTAL * 
* 



* 


qoo 


=152 


0000 




* 
* 


Q10 




ZERO 


Qll 


* 
* 


Q20 


= 41 


9283 


021 


* 
* 


Q30 




ZERO 


0.31 


* 
* 


Q40 


= 4 


8357 


Q41 



ZERO 

ZERO 022 = 0.1730 

ZERO Q32 = ZERO Q33 = ZERO 
ZERO Q42 = 0.0074 Q43 = ZERO 

* 

************************************************************** 

************************************************************** 
* 

* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS 
* 

* R_RMS = 5.5449 X_RMS = 2.3841 Y_RMS = 2.3703 Z. 

* R_GEO = 7.1584 X_GEO = 5.3310 Y_GEO = 5.3001 Z. 
* 

************************************************************** 
************************************************************ 

* BOHR DEFORMATION PARAMETERS (FIRST-ORDER APPROXIMATION) 
* 

************************************************************** 
* 

* BIO = ZERO Bll = ZERO 

* 

* B20 = 0.7112 B21 = ZERO B22 = 0.0029 

* 

* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO 
* 

* B40 = 0.5075 B41 = ZERO B42 = 0.0011 B43 = ZERO 
* 



Q44 =-9.2E-04 * 
* 

***************** 

***************** 
* 

TOTAL * 

* 

RMS = 4.4095 * 
* 

GEO = 9.8599 * 
* 

***************** 
*************** 

TOTAL * 
* 

***************** 



B44 =-1.4E-04 1 

H 
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* ANGULAR MOMENTA AND THE FIRST MOMENTS OF INERTIA FOR OMEGA = 0.500000 MEV * 

* * 
**************************************************************** 



NEUTRONS 
PROTONS 



ORBITAL 



27.72295 
19.01301 



SPINS 
INTRINSIC 



TOTAL 



0.90622 
0.65769 



28.62917 
19.67069 



ORBITAL 



55.44591 
38.02601 



__J_(1)__ 

INTRINSIC 

1.81243 
1.31537 



TOTAL 



57.25834 
39.34138 



* TOTAL 46.73596 1.56390 48.29986 93.47192 3.12780 96.59972 * 

* ^ ^ ^ * 



NEUTRON CONFIGURATIONS 



* 
* 




P 


S 


12 13 14 15 16 17 18 19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


32 


* 
* 


* 
* 
* 


CONF: 


+ 


+ 


11111111 


1 


1 


1 
































* 
* 
* 


* 
* 


VACC: 


+ 


+ 


11111111 


1 


1 


1 
































* 
* 


* 


CONF: 


+ 




11111111 


1 


1 


1 
































* 


* 


VACC: 


+ 




11111111 


1 


1 


1 
































* 


* 
* 


CONF: 




+ 


11111111 


1 


1 



































* 
* 


* 
* 


VACC: 




+ 


11111111 


1 


1 



































* 
* 


* 


CONF: 






11111111 


1 


1 



































* 


* 


VACC: 






11111111 


1 


1 



































* 



******************************************************************************* 
* * 

PROTON CONFIGURATIONS 



* 
* 




P 


S 


7 


8 


9 


10 


11 


12 


13 14 15 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


* 
* 


* 
* 
* 


CONF: 


+ 


+ 


1 


1 


1 


1 


1 


1 


111 


1 



































* 
* 
* 


* 
* 


VACC: 


+ 


+ 


1 


1 


1 


1 


1 


1 


111 


1 



































* 
* 


* 


CONF: 


+ 




1 


1 


1 


1 


1 


1 


111 


1 



































* 


* 
* 


VACC: 


+ 




1 


1 


1 


1 


1 


1 


111 


1 



































* 
* 


* 


CONF: 




+ 


1 


1 


1 


1 


1 


1 


111 


1 


1 
































* 


* 
* 


VACC: 




+ 


1 


1 


1 


1 


1 


1 


111 


1 


1 
































* 
* 


* 


CONF: 






1 


1 


1 


1 


1 


1 


111 


1 


1 
































* 


* 


VACC: 






1 


1 


1 


1 


1 


1 


111 


1 


1 
































* 



******************************************************************************* 

* * 

* ENERGIES (MEV) * 

* * 
******************************************************************************* 



* 


KINETIC: 


(NEU)= 


1651 


194254 


(PRO)= 


1107 


180948 


(TOT)= 


2758 


375202 


* 


SUM EPS: 


(NEU)= 


-2033 


727026 


(PRO)= 


-1132 


746046 


(TOT)= 


-3166 


473072 


* 
* 


PAIRING : 


(NEU)= 





000000 


(PRO)= 





000000 


(TOT)= 





000000 


* 
* 


COULOMB : 


(DIR)= 


564 


274457 


(EXC)= 


-25 


350821 


(TOT)= 


538 


923637 


* 


CONSTR. 


(MULT)= 





000051 


SLOPE= 





001433 


CORR.= 


-0 


030046 


* 


CONSTR. 


(SPIN)= 


-24 


149930 


SLOPE= 





500000 


CORR.= 


-12 


074965 


* 
* 


REARRANGEMENT ENERGY FROM THE 


SKYRME DENSITY - 


-DEPENDENT 


TERMS= 


1008 


649255 


* 
* 


ROUTHIAN 


(TOTAL 


ENERGY PLUS 


MULTIPOLE 


AND SPIN CONSTRAINTS) = 


-1233 


193474 


* 


SPIN-ORB 


(EVE)= 


-80 


525431 


(ODD)= 


-0 


441965 


(TOT)= 


-80 


967396 


* 
* 


SKYRME: 


(EVE)= 


-4503 


848854 


(ODD)= 


-2 


493580 


(TOT)= 


-4506 


342434 


* 
* 


TOTAL : 


(STAB)= 





000143 


(SP) = 


-1209 


043452 


(FUN)= 


-1209 


043595 



t REAL-CLOCK EXECUTION TIMES IN SUBROUTINES t 

* INCLUDING THEIR COMPLETE DOWN-CALLING TREES * 

* * 
******************************************************************************* 



1023 
1 
3 

261 
94 
24 



=> HFODD 
=> NILSON 
=> AVPARI 
=> DENSHF 
=> INTCOU 
=> INTMAS 



=> GEOMFC 
3 => INTKIN 
96 => NILABS 
64 => DENMAT 
16 => SKFILD 
12 => INTCEN 



21 => POWALL 

180 => DIASIG 

64 => AVIMRE 

37 => MOMETS 

7 => RECORD 

45 => INTSOR 



=> RECOUL 

174 => DIAMAT 

37 => INTMUL 

50 => BEGINT 

153 => INTEGH 

8 => ANGYSP 
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